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Abstract. We use the radiation transfer simulation of Xilouris et al. (1999) to constrain the quantity of dust in 
three nearby spiral galaxies (NGC 4013, NGC 5907 and NGC 4565). The predicted visual optical depth from the 
model is compared with the thermal continuum radiation detected from NGC 4013 and NGC 5907 at 850 [im and 
from NGC 4565 at 1.2 mm. The former is based on SCUBA images of NGC 4013 and NGC 5907, reduced and 
presented for the first time in this work. The comparison of visual optical depth and 850 fim (1.2 mm) emission 
yields the emissivity of dust grains in the submillimeter (millimeter) waveband. We infer a value of 1.2 x 10 4 for 
the emissivity at 850 /im which is a factor 4 higher than the benchmark, semi-empirical model of Draine & Lee 
(1984). At 1.2 mm our values are a factor 1.5 higher than this model. Our estimates are more closely aligned with 
recent measurements carried out in the laboratory on amorphous carbon and silicate particulates. A comparison 
between the distribution of 850 /im (1.2 mm) surface brightness and the intensity levels in the 12 CO(1-0) and 21 
cm lines underlines the spatial association between dust detected in the submillimeter/millimeter waveband and 
molecular gas clouds. We suggest that the relatively high emissivity values that we derive may be attributable to 
amorphous, fluffy grains situated in denser gas environments. 

Key words. ISM: dust, extinction, Galaxies: individual: NGC 4013, NGC 5907, NGC 4565, Galaxies: ISM, Galaxies: 
spiral, Infrared: galaxies, Submillimeter 



1. Introduction 

The continuing improvement in Far-Infrared (FIR) instru- 
mentation has permitted a more complete spectral sam- 
pling of the thermal continuum radiation emitted by dust 
grains situated in external spiral disks. The result has been 
that the FIR optical depth of some nearby galaxies is now 
constrained to within a factor 2-3 (Dupac et al. 2003). 
Ironically, the grain parameter allowing the FIR optical 
depth to be converted into dust column density and ulti- 
mately dust mass - the emissivity Q - still remains un- 
certain by nearly an order of magnitude (Hughes et al. 
1993). Q indicates the efficiency with which dust grains 
of a particular temperature emit FIR thermal radiation 
(sometimes this property is expressed as the mass absorp- 
tion coefficient k; both Q and k will be defined formally 
in Sect. 2). Due to a dearth of direct measurements, a 
large proportion of submillimeter (submm) and millimeter 
(mm) astronomers rely on a single estimate of emissivity 
at 125 fim obtained 20 years ago for a single Galactic re- 
flection nebula (Hildebrand 1983). Given, as we shall see, 
that emissivity varies both strongly with wavelength (A) 
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and environment, Q can be predicted to within little bet- 
ter than an order of magnitude at wavelengths close to 1 
mm. Another work frequently cited is that of Draine & Lee 
(1984). These authors propose emissivities on the basis of 
laboratory experiments for A < 60 /jm and use primar- 
ily a model based on solid state theory to extrapolate to 
longer wavelengths. The 100 /im-brightness per H-atom, 
predicted by Draine & Lee (1984) for high-lattitude dust 
clouds, agrees within 30% of the measurements carried out 
by the Infrared Astronomical Satellite (IRAS). However, 
as the authors themselves acknowledge (p. 107 of Draine 
& Lee 1984), the emissivity proposed for A > 300 fj,m 
may be too low by a factor of 3-4 to be consistent with 
astronomical observations. 

In this paper, we carry out a direct measurement of 
the submm/mm emissivity of dust grains situated in three 
nearby, quiescent spiral disks - NGC 4013, NGC 5907 and 
NGC 4565. The technique, explained in detail in Sect. 2, 
compares two tracers of what is believed to be the same 
dust grain population in order to infer the emissivity. 
The tracers in question are submm/mm surface bright- 
ness (Sect. 3) and visual optical depth (Sect. 4). We have 
already undertaken a similar analysis for the nearby spiral 
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NGC 891 and, in this case, our results imply an 850 \im 
emissivity which is 3 times higher than the Draine & Lee 
(1984) value (Alton et al. 2000). In Sect. 7, we compare 
our results with estimates of emissivity in the literature 
for grains situated in various astrophysical environments. 
Apart from estimating the submm/mm emissivity we ex- 
plore how the dust traced in thermal continuum emission 
relates to the various gas phases (H 2 , HI) within the disk 
(Sect. 8). 

2. Technique for determining the emissivity 

Alton et al. (2000) have already discussed the method in 
some detail but we summarise here the basic precepts. 
Classical dust grains of radius 0.1 [im will tend to reach 
an equilibrium temperature T when immersed in a stel- 
lar radiation field. Under such circumstances, heating due 
to the absorption of optical, ultraviolet and near-infrared 
photons is exactly balanced by cooling due to emission of 
mid and far-infrared radiation. The efficiency with which 
grains emit submm/mm radiation depends on their com- 
position and structure and can be expressed by the emis- 
sivity Q, a dimensionless quantity which indicates how the 
flux density (F) recorded at wavelength A compares with 
that emitted by a blackbody: 



F{\) = jpQ(\)B(\,T) 



(1) 



Here D is the distance to a dust cloud containing n dust 
grains of geometrical cross-section a. B is the planck func- 
tion for a blackbody of temperature T. 1 

In the past, a major obstacle has been constrain- 
ing B(X,T) which, near the peak of the FIR spectrum 
(A ~ 200 fxm for spiral galaxies), is highly sensitive to the 
adopted grain temperature T. The temperatures we shall 
use, however, are derived from fitting the Raylcigh- Jeans 
tail of spectrum and in this regime B ccT approximately. 
We shall argue that moderate errors in T due to the fitting 
process (~30%) are not the major source of uncertainty in 
deriving Q. Inspecting Eq. (1), we recognise that, even if 
F and B are known with relative certainty we are required 
to disentangle no from Q. We achieve this by substituting 
the visual optical depth defined as follows: 



T V = NoQ{V) 



(2) 



where Q(V) is the extinction efficiency in the V-band 
(A = 0.55 /im) and N is the number of grains per unit 
area on the sky. At the same time we substitute the sur- 
face brightness /(A) for the flux density F(X) so that 
the distance to the dust cloud D is eliminated. This re- 
formulation of Eq. (1) yields: 



TV 



(3) 



Q(V) 2.2 x 10~ 18 
W)~ QW B(\T ) 

1 A related quantity is the dust absorption coefficient k 
which, for grains of radius a and material density p, is related 
to Q as follows: k = J^2-. 



Table 1. Best-fit parameters derived from RT-modelling of 
NGC 4013, NGC 5907 and NGC 4565 in the V-band. t v is the 
visual optical depth when the galaxy is viewed face-on and 9 
is the inferred inclination of the object, hd and Zd are, respec- 
tively, the exponential scalelength and exponential scaleheight 
of the dust assuming the distance listed in the second column. 
Distances to NGC 4013, NGC 4565 and NGC 5907 are taken, 
respectively, from Bottema (1995), Rupen (1991) and van der 
Kruit & Searle (1982). 



Galaxy Distance 
(Mpc) 



77 



hd 



V ,h d z d v 

(kpc) (kpc) (deg. 



NGC 4013 
NGC 4565 
NGC 5907 



12 
10 
11 



0.67 
0.63 
0.49 



2.45 
9.7 
5.29 



0.13 
0.21 
0.11 



89.7 
88.0 

87.2 



where /(A), for reasons apparent later, is expressed as 
Jy/16" beam. 

In general, it is the ratio qtjt which is of immedi- 
ate interest in comparing our results with other stud- 
ies. However, if we are to quantify the number of FIR- 
cmitting grains in Eq. (1), and thereby ultimately deter- 
mine the dust mass of the system, we should derive the ab- 
solute emissivity Q(X). To do this we must specify Q(V), 
B(X,T), /(A) and t v in Eq. (3). The first of these quanti- 
ties, Q(V), is believed to lie, with some certainty between 
1 and 2 (Whittct 1992, (p. 60); Alton 1996) and we assume 
a value of 1.5 hereafter. The temperature associated with 
the blackbody B(X,T) follows from the shape of the FIR 
spectrum. The submm/mm surface brightness, /(A), and 
the optical depth ry are the two quantities which we shall 
bring together in this paper in order to determine Q(X). 
We start by discussing the submm/mm surface brightness 
and in the subsequent section we show how tv is derived. 

3. Submillimeter/Millimeter Data 

Our submm/mm data originate from the following 
sources. The archive of the James Clerk Maxwell 
Telescope (JCMT) yielded 450 and 850 /im exposures 
of both NGC 4013 and NGC 5907. These raw images 
had been taken with the Submillimcter Common User 
Bolometer Array (SCUBA) and have never been pub- 
lished previously. A reduced and fully-cleaned image of 
NGC 5907 at 1.2 mm, taken with the bolometer array at 
IRAM (Institut de Radioastronomie millimetrique) , was 
provided by Dumke et al. (1997). For NGC 4565, we ex- 
tracted the IRAM major axis profile at 1.2 mm which 
appears in Neininger et al. (1996). 

SCUBA is mounted at the Nasmyth focus of JCMT 
and provides simultaneous imaging at 450 and 850 \xm 
for a region of sky 2.3' in diameter (Holland et al. 1999). 
The shortwave array (450 fxm) consists of 91 bolometers 
(HPBW = 7.5") whilst the longwave array (850 /im) is 
composed of 37 elements (HPBW = 14.7"). In order to 
provide fully sampled images the secondary mirror moves 
in a 64-step jiggle pattern with the integration time last- 
ing 1 sec at each position. At the same time the secondary 
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Fig. 1. NGC 4013 (top) and NGC 5907 (bottom) observed in 
the 850 ^m continuum with SCUBA. Both the background im- 
age and the contours refer to emission at 850 /mi. For NGC 
4013 the contour levels are at 10.0, 10.4, 11.4, 13.3 15.9, 19.2, 
23.3, 28.1, 33.7 and 40.0 mjy/16" beam while for NGC 5907 
the contour levels are at 20.0, 20.8, 23.4, 27.8, 33.8, 41.6, 51.1, 
62.3, 75.3 and 90.0 mjy/16" beam. In both cases the first con- 
tour level is at the 3 a level and the beam size is 16"(FWHM). 



mirror chops at 7 Hz so as to remove the reference (sky) 
emission. After the 16 steps of jiggle pattern, the telescope 
nods in order to allow for slowly varying sky gradients. For 
both NGC 4013 and NGC 5907 the chop throw was 120" 
perpendicular to the major axis. Hourly 'sky-dip' mea- 
surements allow the sky opacity to be monitored during 
the course of the eight-hour observing shift. 

The dedicated SCUBA software package, SURF 
(Jenness et al. 1997), was used to clean, flat- field and cal- 




Fig. 2. Radiation transfer modelling of NGC 4565 in the V- 
band. The real galaxy, as observed by Howk & Savage (1999), 
is shown at the top whilst the corresponding image generated 
by our radiation transfer simulation is situated at the bottom. 



ibrate the raw images of NGC 4013 and NGC 5907 ac- 
cording to atmospheric attenuation. SURF also facilitates 
a reduction in the final image noise by compensating for 
spatially-correlated sky emission across the field of view. 
The pointing stability of the observations was better than 
5" and the typical 850 /jm zenith opacity was 0.16 for 
NGC 4013 and 0.34 for NGC 5907. 

To complement the SCUBA and IRAM images we ob- 
tained resolution-enhanced IRAS images of NGC 4013, 
NGC 5907 and NGC 4565 at 60 and 100 \xm (HiRes im- 
ages made available by IPAC). The mean spatial resolu- 
tion of these data was ~ 1' at 60 /im and ~ 1.5' at 100 /im. 
The SCUBA 850 /xm images for NGC 4013 and NGC 5907 
are shown in Fig. 1. The corresponding instrumental beam 
is ~ 16" (FWHM). The SCUBA 450 /im data are much 
noisier than the longwave measurements and are not pre- 
sented here. For NGC 4013, we were still in a position to 
obtain a global flux density at 450 /im using the precepts 
outlined in Sect. 5. 

For both NGC 4013 and NGC 5907 the 850 /mi mor- 
phology is fairly clumped. NGC 4013 possesses an asym- 
metric distribution with respect to the nucleus with the 
north-east side of the disk distinctly brighter. An enhance- 
ment in 850 /im emission from one side of the disk was a 
feature noted in the morphology of NGC 891 (Alton et 
al. 1998). Such 'hotspots' cannot be attibuted to viewing 
effects (e.g. a spiral arm located on the near-side) since 
the emission in the FIR is expected to be optically-thin. 
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Table 2. Global flux densities, as defined in Sect. 5 for NGC 
4013, NGC 5907 and NGC 4565. For each galaxy, the flux 
density F\ is given in Jy at wavelength A (in microns). The 
uncertainty in the flux density is ± 15%. 



Galaxy 



F 60 

(Jy) 



Fioo 

(Jy) 



F450 
(Jy) 



F850 

(Jy) 



F1200 
(Jy) 



NGC 4013 
NGC 4565 
NGC 5907 



7.0 
7.1 
16 



23 
29 
56 



4.1 



0.61 



1.6 



0.55 
0.54 



Table 3. Dust components derived from SED fitting of NGC 
4013, NGC 5907 and NGC 4565 in the FIR/submm/mm 
regime. We decompose the spectrum into 2 components; warm 
dust of temperature T w and cold dust of temperature T c . The 
ratio of cold dust mass to warm dust mass is given as -£™- . An 
index of 1.5 has been assumed for the wavelength dependency 
of the emissivity (Q ex A -1 ' 5 ). 



Galaxy 



c 



ryW rp\ 

(K) (K) 



M"' 



NGC 4013 39 
NGC 4565 43 
NGC 5907 28 



23 49 
20 99 
13 5.7 



4. Derivation of visual optical depth 

Maps of visual optical depth for NGC 4013, NGC 5907 
and NGC 4565 were derived using the Radiation Transfer 
(RT) model of Xilouris et al. (1999). This numerical model 
simulates the absorption and scattering of optical and 
near-infrared photons by dust grains situated in a stel- 
lar disk using the Henyey-Greenstein phase function and 
the Galactic albedo (Lillie & Witt 1976). Dust and stars 
are assumed to be mixed and smoothly distributed. Both 
components behave independently, falling off as exponen- 
tial distributions in both z-height and radius. The bulge is 
treated as a dust-free component containing stars decreas- 
ing as Ri . The scalelength, scaleheight and optical depth 
of the exponential dust disk are adjusted until there is a 
match between predicted surface brightness and observed 
surface brightness over the entire edge-on disk. 

Our method is best suited to galaxies lying within 5° 
of edge-on (hence the current sample of highly-inclined 
galaxies). This flat orientation allows us to: (1) approx- 
imate, within a limited extent, a clumpy dust disk by a 
smooth exponential distribution and (2) use the unatten- 
uated surface brightness at large z-heights as a boundary 
condition for fitting the dust distribution (see Fig. 1 in 
Xilouris et al. 1997). Although the model possesses certain 
limitations (to be discussed in Sect. 7.4), the opacities in- 
ferred independently in several wavebands (B, V, R, I, J, 
K) are consistent with a Milky Way extinction law within 
a few percent giving us a fairly high degree of confidence 
in the simulation (see Fig. 12 of Xilouris et al. 1999). 

The RT simulations for NGC 4013 and NGC 5907 have 
already been presented in Xilouris et al. (1999) and we 
adopted the corresponding dust parameters without fur- 
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Fig. 3. Flux densities for NGC 4013, NGC 5907 and NGC 
4565 as observed by IRAS (60 and 100 fim), SCUBA (450 
and 850 (im) and IRAM (1.2 mm). In each case, the thermal 
spectrum has been fitted by the sum of a warm and cold dust 
component (dashed and dotted lines respectively). The total 
SED is given by the solid curve. The proportion of dust in each 
of the warm and cold components, as well as their associated 
temperatures, are given in Table 3. 
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ther treatment. The simulation for NGC 4565 was new 
and required at least one optical/NIR image for the fit- 
ting procedure. To this end, a reduced and cleaned V- 
band image of NGC 4565 was retrieved from the NED 
server (originally taken by Howk &; Savage (1999) using 
the WIYN telescope). Fig. 2 illustrates both the observed 
and simulated images for this object. Table 1 lists the 
model parameters, including the visual optical depth, for 
all three galaxies in the sample. As indicated in this table, 
we adopt distances of 12, 10 and 11 Mpc to NGC 4013, 
NGC 4565 and NGC 5907, respectively. 

5. Grain temperature from SED fitting 

The spectral energy distribution (SED) for NGC 4013 
and NGC 5907 was derived in the following manner. The 
SCUBA, IRAM and HiRes 60 /Ltm images were smoothed 
to a common 'worst' spatial resolution namely that of the 
100 /im IRAS image (~ 1.5' FWHM). To avoid aperture 
corrections (Chini et al 1995; Hughes et al 1997) an image 
area was defined over which the object had been detected 
in all FIR/submm/mm filters. The emission within this 
area was used as the global flux density. The arcminute 
resolution of the IRAS data is too coarse to monitor the 
change in grain temperature with galactocentric radius 
and, therefore, we sought only to fit global mean dust 
temperatures to the thermal spectrum of each object. For 
NGC 4565 we integrated the major axis profile at 1.2 mm 
appearing in Neininger et al. (1996) and measured the 60 
and 100 (im flux densities over the same area of our HiRes 
images. The global flux densities we derive for all three 
objects are shown in Table 2 and the corresponding SEDs 
are depicted in Fig. 3. 

To determine the grain temperature we assume that 
the flux densities recorded in Fig. 3 arise purely from ther- 
mal continuum emission. Contamination of the 450 and 
850 ^m filters by the respective 12 CO(6-5) and 12 CO(3- 
2) line transitions was neglected in the first instance. 
Similarly, no account was taken at this stage of non- 
equilibrium emission from small grains. Thus, we decom- 
posed each SED into two dust components (warm and 
cold) assuming that each component emits as a modified 
blackbody with (3 = 1.5 [F(X,T) ex B(\,T)\-P}. Our rea- 
sons for chosing this relatively low value of [3 will become 
evident at a later stage in the paper (Sect. 7.2). The warm 
and cold grain temperatures derived from the fitting pro- 
cedure are summarized in Table 3. We emphasize that our 
aim here is to fit the shape of the SED. We are not seek- 
ing to derive absolute dust masses which would require 
adopting a priori a FIR emissivity. We derive 3 parame- 
ters (cold dust temperature, warm dust temperature and 
ratio of cold to warm dust mass) from either 3 (NGC 4565) 
or 4 data points (NGC 4013, NGC 5907). 

From Table 3, we surmise that 85-99% of the dust de- 
tected is cold possessing a temperature of 13-23 K. Similar 
values are recorded for diffuse Galactic dust mapped with 
COBE (Reach et al. 1995) and from authors carrying 
out more complete spectral mapping of external galaxies 
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Fig. 4. Major axis profiles of NGC 4013, NGC 5907 and NGC 
4565 in 850 /im and 1.2 mm surface brightness (solid line, left 
axis) and visual optical depth (dashed line, right axis). The 
latter has been inferred from radiation transfer modelling and 
the map in visual optical depth has been smoothed to the cor- 
responding spatial resolution of the submm/mm data before 
profiling (16" and 20" at 850 /xm and 1.2 mm respectively). 
The major axis is defined such that positive values correspond 
to north-east, south-east and north-west for NGC 4013, NGC 
5907 and NGC 4565 respectively. 
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(e.g. Dupac et al. 2003). A discussion of the uncertain- 
ties in the grain temperature and its implications for the 
inferred emissivity is reserved to Sect. 7.2. 



6. Deriving the emissivity Q 

In order to compare our two chosen tracers of galactic 
dust - visual optical depth and submm/mm flux density 
- on equal terms we are compelled to smooth the map of 
Ty from the RT model to the same spatial resolution as 
the SCUBA/IRAM data from which we wish to determine 
the emissivity. For NGC 4013 and NGC 5907 our analy- 
sis was for the 850 /im waveband whilst for NGC 4565 
the comparison was conducted at 1.2 mm. We emphasize 
that the 850 /im and 1.2 mm filters lie on the Rayleigh- 
Jeans tail of the SED which renders our comparison less 
susceptible to uncertainties in the adopted grain tempera- 
ture. After smoothing the map in V-band optical depth to 
the required spatial resolution we profiled along the ma- 
jor axis of the galaxy in order to derive the mean ratio 
-I7JY [required for Eq. (3)] . These profiles are illustrated in 
Fig. 4. Profiles were also generated perpendicular to the 
major axis and these indicate that visual extinction and 
850 /im surface brightness decrease in a similar manner 
for increasing z-height (distance above the midplane) . We 
do not attach too much significance to the correlation in 
z-height because our spatial resolution (16") is too coarse 
with respect to the expected exponential scaleheight of 
galactic dust (0.25 kpc or ~ 5" for a galaxy 10 Mpc away). 

The profiles in Fig. 4 indicate that the visual optical 
depth falls off radially in a manner similar to that of the 
submm/mm surface brightness. We have, therefore, some 
confidence that the two quantities are tracing the same 
component of galactic dust. Spikes in the submm/mm 
emission are evident which cannot be accounted for by 
the smooth dust distribution adopted by the RT model. 
We emphasize that such localized peaks, at least for NGC 
4013, are real within the noise and, as such, might either 
correspond to concentrations of molecular gas or warmer 
thermal emission from grains situated within star-forming 
regions. A discussion of how the submm/mm emission re- 
lates to the various gas phases of the ISM is deferred to 
Sect. 8. The percentage of 850 /im emission contained in 
the spikes of Fig. 4 can be gauged by fitting a smoothly- 
declining baseline to the profile and totalling the emission 
above this level. The method is somewhat crude but yields 
a percentage of ^20% for the locally-enhanced 850 /im flux 
density in NGC 4013. We conclude that both the visual 
extinction and submm thermal emission appear to trace 
the same component of galactic dust but a significant frac- 
tion (~ i) of the submm emission cannot be accounted 
for by our simplified model of visual extinction. 

From the profiles in Fig. 4 we obtain mean values of 



and 1.2 mm in Jy/16" beam and Jy/20" beam respec- 
tively. In Sect. 5 we have decomposed the SED of each 
object into a warm and a cold dust component in or- 
der to determine the grain temperature for the black- 
body function. We now intend to use Eq. (3) to deter- 
mine Q(850 pm) and Q(1.2 mm) but, before doing so, 
we subtract off the minor contribution (< 20%) of the 
warm dust component to /(850/im) and /(1.2 mm). This 
allows us to substitute a single grain temperature into 
Eq. (3) namely that of the cold dust component listed in 
Table 3. This simplified approach yielded a ratio — ^ ' 



51 ± 2 and 50 ± 2 for the ratio 



in NGC 4013 



/(850^m) 

and NGC 5907, respectively, and a value of 330 ± 40 
for the ratio /(1 ^ m) in NGC 4565. Here, /(850 pm) 
and /(1.2 mm) refer to the surface brightness at 850/im 



Q(850^m) 

of 14000 and 12000 for NGC 4013 and NGC 5907, re- 
spectively, and a ratio f{ ^ m} = 66000 for NGC 4565. 
Assuming Q(V) ~ 1.5 (Sect. 2), we infer an emissivity of 
~ 1.2 x 10~ 4 and ~ 2.3 x 10~ 5 at 850 /im and 1.2 mm, re- 
spectively. For our adopted grain parameters (a = 0.1 /im 
and p = 3000 kgm -3 ) this equates to 3.0 and 0.57 cm 2 g _1 
for k(850 jim) and k(1.2 mm) respectively. 



7. Discussion on the emissivity (Q) 

In Table 4 we compare our results with existing estimates 
of submm/mm emissivity for a diversity of astrophysical 
media (diffuse HI clouds, reflection nebulae etc). The cor- 
responding dust absorption coefficient k is given assuming 
a classical grain size of 0.1 /im, Q(V) = 1.5 and a material 
density typical of amorphous silicates p = 3 x 10 3 kgm~ 3 
(Mcnnella et al. 1998). The same data are presented 
graphically in Fig. 5 with identical markers denoting simi- 
lar grain environments. Most of the points in Fig. 5 rely, as 
we do in this paper, on observed optical/NIR extinction to 
quantify the amount of dust emitting in the submm/mm 
regime. Notable exceptions are the semi-empirical model 
of Draine & Lee (1984) and the recent analysis of James 
et al. (2002). The former, already alluded to in Sect. 1, 
constitutes a benchmark study which is frequently cited 
in the literature (see also Li & Draine (2001), a recent 
update of the model which proposes only modestly differ- 
ing emissivities in the FIR). James et al. (2002) derives 
the emissivity by quantifying the amount of dust on the 
basis of local gas density and local metal abundance (a 
constant fraction of metals is assumed to be tied up in 
dust grains). In a similar manner, estimates by Dumke et 
al. (1997) and Neininger et al. (1996) for NGC 5907 and 
NGC 4565, respectively, are based on the column density 
of atomic hydrogen detected in the outer regions of the 
galaxy. A solar gas-to-dust ratio is employed in order to 
quantify the amount of dust present and the level of 1.2 
mm emission is compared with this quantity of dust in 
order to infer the emissivity. 

The values derived in this paper, plus that of NGC 891 
based on the same technique, indicate a low ratio Q?p IR \ 
and therefore a high efficiency for submm/mm thermal 
emission compared with, for example, COBE measure- 
ments of high-lattitude dust in the Milky Way (difference 
of factor 4 at 850/im). Indeed, our estimates are more 
aligned with denser gas environments, where grains are 
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Table 4. Values for submillimeter and millimeter emissivity, Q(A), collated from the literature. Estimates are given as q, x < where 
the extinction efficiency in the V-band [Q(V)] has been taken, where necessary, to be f .5 (Whittet 1992). The absorption mass 
coefficient k(X) is calculated from Q(X) assuming a classical grain radius of 0.1 ^irn and a material density of p = 3 x 10 3 kgm~ 3 . 
For each measurement the column 'Medium' denotes the type of astrophysical environment where the corresponding grains are 
believed to be situated. 



A 
(H 



~qwr 

Q(A) 



k(A) 

(cm'g- 1 ) 



Medium 



Reference 



Comment 



100 


760 


49 


diffuse HI (Milky Way) 


Bianchi et al. (1999) 


125 


2000 


19 


reflection nebula 


Hildebrand (1983) 


240 


3900 


9.7 


diffuse HI (Milky Way) 


Boulanger et al. (1996) 


250 


2500 


15 


reflection nebula 


Casey (1991) 



250 



4400 



8.5 



diffuse HI (model) 



400 


1900 


20 


late-type star 


736 


34000 


1.1 


diffuse HI (Milky Way) 


850 


9100 


4.1 


cold core 


850 


30000 


1.3 


cold core 


850 


54000 


0.69 


diffuse HI ? 


850 


16000 


2.4 


H 2 (NGC 891) 


850 


14000 


3.2 


H 2 (NGC 4013) 


850 


12000 


2.7 


H 2 (NGC 5907) 


850 


25000 


1.5 


cold core 


1000 


18000 


2.1 


diffuse HI (model) 


1000 


180000 


0.22 


laboratory (Si crystalline) 


1000 


40000 


0.95 


laboratory (Si amorphous 


1000 


1600 


24 


laboratory (C amorphous) 


1200 


120000 


0.315 


diffuse HI (NGC 4565) 


1200 


110000 


0.35 


diffuse HI (NGC 5907) 


1200 


110000 


0.34 


cold core 


1200 


66000 


0.57 


H 2 (NGC 4565) 



Draine & Lee (1984) 

Sopka (1985) 
Boulanger et al. (1996) 

Rengarajan (1984) 
Kramer et al. (2003) 
James et al. (2002) 
Alton et al. (2000) 
this work 
this work 

Bianchi et al. (2003) 
Mathis & Whiffen (1989) 
Agladze et al. (1996); 
Mennella et al. (1996) 
Agladze et al. (1996); 
Mennella et al. (1996) 
Mennella et al. (1996) 
Ncininger et al. (1996) 

Dumke et al. (1997) 
Bianchi et al. (2003) 
this work 



Q$Q up to 5000 it (3 = 1 
DIRBE data; values also given 
at 100 and 140 rmi 
5 objects; ^Q up to 30000 
depending on the method 
Uses laboratory data; oscillator 
model for A > 60 /im 

FIRAS data; values also given 
at 346, 490, 535 and 1100 /*m 



Using N H , metallicity & SCUBA data 
RT model of this work & SCUBA data 



SCUBA data 

Mean value of both works 
(150000 and 197000) 
Mean value of both works 
(36100 and 43600) 

Authors claim emissivity 2-4 
higher for molecular regions 
Outer regions of the galactic disk 
IRAM data 



expected to coagulate into amorphous particulates rather 
than remain as well-ordered crystalline cores. Our results 
lie within the range suggested by recent tests on amor- 
phous carbon and silicate particulates conducted in the 
laboratory. Submm/mm emissivity is expected to be en- 
vironment sensitive responding directly to grain temper- 
ature and, indirectly (through coagulation), to ambient 
gas density. Recent observational evidence points to a 
threefold increase in Q(600 /J,m), compared to the diffuse 
ISM, for fluffy, coagulated grains forming in dense dust 
clouds (Stepnik et al. 2003). Experiments in the labora- 
tory on amorphous grains predict a fall of factor 2-6 in 
Q(X ~ 1 mm) as the grain temperature lowers from 300 
K to 24 K. Results presented below indicate a strong asso- 
ciation between 850/mr emission and the 12 CO(1-0) line in 
our sample of galaxies. Gas clouds in the Milky Way that 
emit strongly in the CO-line are characterised by densities 
of ~ 10 3 m~ 3 (Scaquist et al. 2004) and are thus more akin 
to dense clouds rather than the diffuse HI clouds analysed 
in the aforementioned COBE study. James et al. (2002) 
claim to derive the same emissivity regardless of the value 



of X used to convert the 12 CO(1-0) line emission of their 
sample into H2 gas masses. This seems difficult to explain 
unless their sample is dominated by Hi-rich, metal poor 
objects rather than the H2-rich galaxies that make up our 
own sample. 

Before accepting the possibility that the relatively high 
emissivity values we derive are attributable to grains in 
denser H2 clouds, we carry out a critical examination of 
the assumptions and modalities of our technique. These 
can be summarized as follows: (i) the grains responsible for 
optical/NIR extinction are not, as supposed, the origin of 
the submm/mm emission, (ii) our SED fitting is erroneous 
implying an incorrect grain temperature (iii) we overesti- 
mate the submm/mm flux densities (iv) we underestimate 
the visual opacity. We discuss each of these points in turn 
below. 

7.1. Error in the basic assumptions 

Our central assumption is that grains extinguishing opti- 
cal/NIR radiation in spiral disks will also be the main 
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Fig. 5. Emissivity Q(X) of astrophysical dust in the far-infrared, subnim and millimeter wavebands. The plot depicts the values 
collated in Table 4 employing identical markers for similar astrophysical environments. The estimates that we derive for NGC 
4013, NGC 4565, NGC 5907 and NGC 891 (the latter in Alton et al. (2000)) have been attributed to dust residing in molecular 
gas clouds (solid circles). The dashed and dotted lines correspond to a wavelength dependency for the emissivity which varies as 
/?=2.0 and /?=1.5 respectively (Q oc A ). The lines have been set to pass arbitrarily through the frequently-cited model value 
of Draine & Lee (1984) at 250 /urn. 



emitters of submm/mm thermal radiation. Hildcbrand 
(1983) has already shown that the dependencies of opti- 
cal extinction and FIR emission on grain size are broadly 
similar. We make a critical appraisal of this argument in 
the Appendix using the MRN grain-size distribution as 
a test-case (Mathis et al. 1977). To summarize the con- 
clusions of that section, we find that, for the same grain- 
size distribution in all parts of the ISM, grains of roughly 
the same size will dominate both optical extinction and 
submm/mm thermal emission. The enhanced emissivity 
we measure in the previous section, however, may sug- 
gest that our SCUBA images are sensitive to a population 
of grains that are significantly larger than the classical 
0.1/im grains dominating optical extinction. This ques- 
tion of bi-modality in the grain size can only be solved 
expeditiously by a decomposition of galactic submm/mm 
emission into discrete and diffuse sources (as has been 



carried out for the Milky Way, at shorter wavelengths, 
using IRAS). A project of this kind is probably beyond 
the reach of current submm/mm arrays where the field 
of view (2 — 3') and, therefore, ill-suited to a large-scale 
survey of the Galactic plane. Prior subtraction of discrete 
sources from our SCUBA images, in order to eliminate 
larger, non-classical grains, would act to lower the emis- 
sivity we derive for the classical population dominating 
the optical extinction. 

7.2. Uncertainties in SED fitting 

SED fitting is rendering particularly uncertain by the 
value chosen for wavelength-dependency /?. SED fits with 
(3 = 2, rather than the f3 = 1.5 adopted in Sect. 5, 
imply lower grain temperatures and therefore higher 
submm/mm emissivities. For example, we derive 13 K for 
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the cold dust component in NGC 4013 (instead of 23 K 
with (3 = 1.5) and consequently Q /g 50 m) = 6500 (i.e. a 
submm emissivity which is twice as high). In conclusion, 
adopting values of [3 higher than 1.5 would increase the 
divergence of our results from COBE-based estimates of 
emissivity and the Draine & Lee (1984) model. 

Non-equilibrium emission from very small grains and 
PAHs is also expected to contribute to the flux density de- 
tected in the IRAS filters. Consequently we subtract 62% 
and 14% from the flux densities at 60 and 100 ^.m (Desert 
et al. 1990) and re-fit the modified SED to gauge the mag- 
nitude of this effect. The implied change in the emissivity 
is ~ 10%. In summary, since the 850 ^.m filter samples 
emission on the Rayleigh-Jeans tail, the emissivity we in- 
fer at that wavelength is only linearly-proportional to the 
uncertainties in the grain temperatures. Our lattitude in 
fitting T is too small to account for the relatively high 
emissivities we derive. 



7.3. An overestimate of submm /mm emission ? 

The SCUBA filter at 850 ^m and the IRAM contin- 
uum bandpass at 1.2 mm are expected to contain some 
line emission from the respective transitions 12 CO(3-2) 
and 12 CO(1-0). For the former, contamination is expected 
to be as high as 50% for dense gas situated within 1-2 
kiloparsec of the galactic nucleus (Bianchi et al. 2000). 
Similarly, Seaquist et al. (2003) estimate a mean contam- 
ination of 25% for the infrared-luminous galaxies making 
up the SCUBA Local Universe Survey (SLUGS; Dunne et 
al. 2000). For the general disk, however (appropriate to 
the analysis we are carrying out here), the CO gas is ex- 
pected to be less collisionally excited and the correspond- 
ing 12 CO (3-2) contamination is believed to be much closer 
to 10% (Israel et al. 1999; Boettner et al. 2003). Images 
of NGC 5907 in the 12 CO (3-2) emission line (Dumke et 
al in prep.) indicate a contamination of 130mJy to the 
SCUBA continuum filter (i.e. less than 10% of the contin- 
uum level). 

IRAM heterodyne and continuum measurements of 
nearby, quiescent disks imply corrections of ~ 10% are 
necessary to eliminate the 12 CO(1-0) line emission from 
the 1.2 mm continuum filter (e.g. Guclin et al. 1993). The 
contribution of the 12 CO(6-5) line to the SCUBA filter 
at 450 /im is expected to be less than 10% (Seaquist et 
al. 1994). To gauge the impact of non-thermal emission 
on our results we have re-fitted the SEDs in Fig. 3 remov- 
ing, beforehand, both non-equilibrium emission from small 
grains (Sect. 7.2) and a line-contamination of 10%, 20% 
and 10% to the wavebands 450 /xm, 850 /xm and 1.2 mm re- 
spectively. The impact on our infered values of — ^ - 



10 



Q(850/jm) 

is comparatively small (~15%) and cannot account for 
the relatively high submm/mm emissivities which we have 
found. 
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Fig. 6. Dust and neutral gas profiles along the major axis of 
NGC 4013 and NGC 5907. The solid circles denote the distri- 
bution of grains (using multiplicative factors of 500 and 350 to 
clarify the plot). The open circles and open squares trace the 
respective distributions of molecular and atomic gas within the 
disk. The major axis is defined such that positive values cor- 
respond to north-east and south-east for NGC 4013 and NGC 
5907 respectively 

7.4. An underestimate of optical opacity ? 

In a synopsis of disk opacity, Alton et al. (2001) found 
that the optical depth implied by the RT-modcl of Xilouris 
et al. (1999) was generally low compared to other values 
cited in the literature. Many of these other studies, how- 
ever, are based on far-infrared emission (e.g. ISOPHOT 
maps at 200^tm) and assume a priori an emissivity which 
is lower than that derived here (implying necessarily a rel- 
atively high visual optical depth) . This was the conclusion 
of Popescu et al. (2000) and Misiriotis et al. (2001), who, 
using the same RT model as us but employing the lower 
emissivity of Draine & Lee (1984), explain the detected 
submm/mm emission from external galaxies in terms of 
a second, massive dust component distributed 'thinly' in 
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z-height in a similar fashion to the OB stellar population. 
The 'thin' nature of this disk is imperative - increasing 
the optical depth of the 'standard' dust disk we have as- 
sumed here (scale-height of dust half that of the stars) 
would lead to a flagrant discrepancy between the model 
results and the observed photometry. We ask ourselves 
whether this second dust disk, which would increase the 
visual optical depth by a factor of 3, is realistic. IRAS 
and COBE (DIRBE) images of the solar neighborhood 
indicate a dust scaleheight of 0.13 kpc with the possibil- 
ity of cold dust having an exponential scaleheight as great 
0.5 kpc (Davies ct al. 1997). This is larger than the 0.09 
kpc scaleheight characterising the OB stellar population 
within the Galaxy. We admit that we cannot rule out the 
presence of a second, massive dust with our current ob- 
servations and that this must be seriously considered as 
an alternative to the relatively high emissivity we have 
derived. Essentially, with a massive thin disk, the near- 
infrared becomes optically thick and the edge-on surface 
photometry is no longer diagnostic. Moreover, the thick- 
ness of the second dust disk approaches that of the NIR 
spatial resolution. A discriminatory test will be to predict 
the azimuthally-averaged colours of face-on spiral disks 
(e.g. B-K) using (i) an RT-model containing a massive, 
thin dust layer (low emissivity) and (ii) an RT-model with 
only a 'standard' dust layer (high emissivity). A compar- 
ison between model and observation is being undertaken 
by one of us as a thesis project (Dasyra) but no conclu- 
sions from this study are currently available. 

A component of dust in clumps too optically thick 
and compact to be probed by RT simulations is feasi- 
ble. However, as argued by both Alton et al. (2000) and 
Misiriotis & Bianchi (2002), this is unlikely to exceed 50% 
of the mass assumed in the current, smooth dust distri- 
bution. A correction of the magnitude is insufficient to 
explain our results. 

8. The submm as a tracer of gas 

The dust in our RT-models has been assumed to be 
smoothly distributed in order to simplify our calculations. 
The major axis profiles plotted in Fig. 4 indicate, how- 
ever, that a conspicuous fraction of the 850 /xm thermal 
emission detected within our sample is clumped at the 
kpc-level. NGC 4013 emits, for example, much more copi- 
ously in the north-east part of the disk (analogous to the 
enhanced submm emission detected in the north-east disk 
of NGC 891; Alton ct al. 2000). We have argued that a 
possible physical interpretation of the high emissivities de- 
rived in Sect. 7 lies in an association between submm/mm 
emission and concentrations of molecular gas. In this sec- 
tion we examine this hypothesis further and compare the 
submm/mm emission within our sample with the corre- 
sponding distributions of molecular and atomic gas within 
the disk. In this way we hope to learn which medium 
within the disk is harbouring the grains giving rise to the 
submm/mm emission. In the past, tight spatial correla- 
tions have been established between CO line emission and 



the 450 /xm, 850 /zm and 1.2 mm thermal emission de- 
tected from nearby, quiescent disks such as NGC 7331, 
NGC 891 and NGC 6946 (Guelin ct al 1993; Alton ct al. 
2000; Bianchi et al. 2000; Alton ct al. 2001). 

Information regarding the distribution of gas along 
the major axis was obtained from the following sources. 
For NGC 4013, both the atomic and molecular hydrogen 
column densities were extracted from Fig. 5 of Gomez 
de Castro & Garcia-Burillo (1997). For NGC 5907, the 
12 CO(1-0) line intensity and 21 cm emission were taken 
from Fig. 6 of Dumke et al. (1997) using their Fig. 5 to 
calibrate the latter. For an indication of how dust and gas 
are correlated within NGC 4565 we refered to the existing 
study of Ncininger et al. (1996). To convert the 12 CO(1-0) 
emission to a H2 column density we adopted a value of X= 
2 xlO 20 cm~ 2 K km/s (likewise, values of H2 column den- 
sity for NGC 4013 from Gomez de Castro & Garcia-Burillo 
(1997) were adjusted to this value oiX). We recognise that 
there is some debate concerning the exact value of X (es- 
timates vary from about 1.5 to 6 xlO 20 cm~ 2 K km/s for 
large, quiescent spiral disks) but values close to 2 xlO 20 
cm~ 2 K km/s seem appropriate for the general ISM of the 
Milky Way (Maloney 1990). The SCUBA 850 ^m images 
of NGC 4013 and NGC 5907 were smoothed to the same 
spatial resolution as the gas data (20") before comparing 
the distribution of submm/mm emission with the gas. 

Fig. 6 illustrates the major axis distributions of gas and 
dust in NGC 4013 and NGC 5907. The surface mass den- 
sity for the dust follows from the parameter n in Eq. (1) 
assuming the emissivity Q we have derived in Sect. 6 and a 
respective grain radius and material density of 0.1 /im and 
3000 kgrn -3 . The profiles in Fig. 6 reinforce the idea that 
grains emitting in the submm are associated with molecu- 
lar gas, particularly within the central 10 kpc of the disk. 
Moving towards the disk edge, the fall-off in submm emis- 
sion is intermediate between that of the molecular and 
atomic gas for NGC 4013. Ncininger et al. (1996) recog- 
nise a similar trend for NGC 4565 at 1.2 mm and anal- 
yses of the submm and millimeter emission from NGC 
891 present similar radial correlations (Alton et al. 2000; 
Guelin et al. 1993). For NGC 5907, the 850 ^m appears 
to follow the molecular gas at all radii and falls to zero 
when only atomic gas is present. 

The tight correlations we observe between 850 /im 
emission and the 12 CO (1-0) line intensity might lead us to 
the conclusion that the SCUBA longwave filter is contami- 
nated much more severely by 12 CO (3-2) line emission than 
the < 20% we have argued in Sect. 7.3. We emphasize, 
however, that correlations between CO-line emission and 
submm/mm thermal emission have already been noted 
when the continuum filter in use is known to contain only 
minor line contamination. For example, a strong associa- 
tion between 1.2 mm continuum emission and 12 CO (1-0) 
line intensity has already been established in both NGC 
891 (Guelin ct al. 1993) and NGC 5907 (Dumke et al. 
1997). These observations, along with the estimates of CO- 
line contamination produced in Sect. 7.3, persuade us that 
the correlation we find between 850 /jm surface brightness 
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Fig. 7. The gas-to-dust ratio along the major axes of NGC 
4013 and NGC 5907. The CO line has been converted to a H 2 
column density assuming X — 2 x 10 20 cm~ 2 K km/s. The dust 
column density follows from the submm emissivities derived in 
Sect. 6 of this work. 



and 12 CO(1-0) line emission indicate a true physical asso- 
ciation between CO gas and cold dust. 

We use our profiles of gas and dust surface density 
along the major axis to infer the gas-to-dust ratio. This 
is illustrated in Fig. 7 for NGC 4013 and NGC 5907. 
Furthermore, we employ the trend noted by Garnet (1988) 
between absolute B-band magnitude and relative oxygen 
abundance in order to infer the gas-phase metallicity (see 
also Garnett & Shields 1988). Although extinction effects 
are severe in these edge-on disks, our radiation transfer 
model can be used to derive the intrinsic blue luminos- 
ity of NGC 4013 and NGC 5907. The object in question 
can then be located on the O/H vs. Mb plot of Garnet. 
Following Garnet, we adopt Hq=50 kms /Mpc in order 
to determine the distance modulus and hence absolute B- 
band magnitude. In this way, we derive values of -19.8 and 
-20.5 for M B in NGC 4013 and NGC 5907, respectively. 
On the Garnet plot this corresponds to a metallicity of 
z = O.79z (for NGC 4013) and z = z (for NGC 5907) 



i.e. a solar-type metallicity within the 40% uncertainties 
inherent in the estimate. 

The error-weighted gas-to-dust ratios that we infer 
from Fig. 7 are close to 500 and, therefore, a factor 2-3 
higher than the values of 150-300 associated with the so- 
lar neighborhood. We emphasize that our gas-to-dust es- 
timates are sensitive to the chosen value of the conversion 
factor for 12 CO(1-0) to H 2 . Frequently, X ~ 0.5 x 10 20 
cm~ 2 K km/s has been proposed for the collisionally- 
excited climes characterising the central few kpc of spi- 
ral galaxies (e.g. Neininger et al. 1996). Reducing our X 
parameter to this level would bring our gas-to-dust ratio 
much closer to the solar value. 



9. Summary and conclusions 

We have reduced and collated submillimctcr and millime- 
ter imaging data for three nearby edge-on spiral galaxies 
(NGC 4013, NGC 4565 and NGC 5907). We have com- 
pared the distribution and level of thermal continuum 
emission from the dust with the V-band optical depth in- 
ferred from a sophisticated (scattering-l- absorption) radi- 
ation transfer simulation (Xilouris et al. 1999). The anal- 
ysis yields the emissivity of dust at wavelengths of 850 /im 
and 1.2 mm. At 850 /im we infer an emissivity of 1.2 x 10 4 
i.e. a factor 4 higher than both predictions from the Draine 
& Lee (1984) model and measurements made on high- 
lattitude Milky Way dust using COBE. At 1.2 mm, our 
estimate is only a factor 1.5 higher. Our results, particu- 
lary at 850/zm, are consistent with recent measurements 
carried out under laboratory conditions where relatively 
high emissivity values are recorded for amorphous silicate 
and carbon particulates. The submm/mm emission de- 
tected from our objects follows closely the distribution 
of molecular gas, the latter being traced in the 12 CO(1-0) 
line. A physical explanation for the relatively high emissiv- 
ity values we derive might be that dust emitting strongly 
in the submm/mm waveband is situated chiefly in molec- 
ular gas clouds where the elevated density is conducive to 
the formation of amorphous, fluffy grains. Such grains are 
expected to possess emissivity values that are a factor 3 
or so higher than refactory cores circulating in the diffuse 
ISM (Osscnkopf & Hcnning 1994). 

We accept that our radiation transfer simulation could 
'miss' a significant fraction of galactic dust if the grains are 
distributed within a non-standard way within the disk (i.e. 
within a very narrow layer with respect to the mid-plane) . 
A 'litmus test' for both the radiation transfer model, and 
the high emissivities we derive as a consequence, will be 
a prediction of the observed optical and near-infrared 
azimuthally-averaged radial colours of face-on spirals as- 
suming the optical depths we use in this work. For ex- 
ample, an observed B-K colour which is higher than that 
predicted will indicate an underestimate of visual optical 
depth for the current edge-on sample. Work is currently 
being undertaken by one of us (Dasyra, thesis project) to 
resolve this issue. 
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We also accept that a major part of the submm/mm 
thermal emission detected in our SCUBA/IRAM images 
may arise from bigger grains contained within localised, 
compact sources (e.g. pre-stellar cores or circumstcllar 
disks) whilst the optical extinction we measure is almost 
certaintly attibutable to 'classical' grains (characterised 
by a Galactic extinction law). If a bimodality of this kind 
is indeed present in our observations, we would also over- 
estimate the emissivity associated with diffuse interstellar 
dust. 

The gas-to-dust ratios we derive for NGC 4013 and 
NGC 5907 are ~ 500 and therefore a factor 2-3 higher 
than values cited for the solar neighborhood. Our esti- 
mates depend sensitively on the adopted value of X, the 
conversion of 12 CO (1-0) line intensity to H2 column den- 
sity. We have used X = 2 x 10 20 cm~ 2 K km/s in this work 
but, if we were to follow proponents of X — 0.5 x 10 20 cm~ 2 
K km/s, our gas-to-dust estimates would lie close to solar 
values. A crude estimate of metallicity using the Mb-O/H 
relation (Garnet 1988) suggests a metallicity close to solar 
for NGC 4013 and NGC 4565. 

10. Appendix 

Here, we make a critical appraisal of the central supposi- 
tion in this paper, namely that the same grains responsi- 
ble for optical extinction also give rise to the submm/mm 
thermal continuum. The mean radius (a) of grains con- 
tributing to optical extinction T opt can be expressed as 
follows: 



< a >= 



J an(a)T op tda 
J n(a)T opt da 



(4) 



where the integral is weighted by the distribution of grains 
sizes n(a) (e.g. Hildebrand 1983). The linear behaviour 
of the Galactic extinction law within the optical regime 
indicates that Q op t oc j (Whittet 1992). Refering to the 
definition of optical depth [Eq. (2)], Eq. (4) then becomes: 



< a >= 



J an(a)a 3 da 
J n(a)a 3 da 



(5) 



We use the MRN size distribution as a test case where 
n(a) oc a~ 3,5 over the range 0.01 [im > a > 0.25 \xm 
(Mathis et al. 1977). This yields a weighted grain radius 
of 0.1 /im corresponding to the frequently-cited 'classical' 
grain size. We emphasize however that, within the limits 
set out below, the exact size distribution is usually not im- 
portant because emission in the submm/mm can also be 
shown to behave as Eq. (5). For wavelengths much larger 
than the grain size, extinction is dominated by alsorption 
and the Kramcrs-Kronig relations yields the approxima- 
tion Q oc a in the submm/mm regime (Whittet 1992). It 
is the same parameter Q which determines the efficiency 
with which grains emit submm radiation. Thus, for a sin- 
gle grain population characterised by an MRN size distri- 
bution, we expect grains of size 0.1/im to dominate both 
optical extinction and submm/mm thermal emission. 



The arguments advanced above are less robust if grains 
of different sizes are believed to dominate different parts 
of the ISM. Thus, we can envisage a scenario whereby 
optical extinction arises from a widely-distributed popu- 
lation of classical dust grains whereas a significant fraction 
of submm/mm thermal emission emanates from modified, 
larger grains situated in dense, cold gas regions. There is 
some evidence that a bimodality of this kind has been de- 
tected in the present study. The optical/NIR extinction 
law within our sample of galaxies is very similar to that of 
the Milky Way (Sect. 4) indicating an absolute-to-selective 
extinction ratio of Ry — 3. In contrast, the emissivity 
we measure in the submm/mm is relatively high pointing 
to grains of cither an amorphous structure or larger size. 
There is some evidence that dark dust clouds, possessing 
elevated submm/mm emissivities, are also characterised 
by higher total-to-selective extinction ratios (Ry = 4 — 6) 
and hence larger grain sizes (Kandori et al. 2003). 

Bimodality would not normally be expected to be a se- 
vere problem in this study if the two grain populations in 
question are well mixed and uniformly distributed with re- 
spect to the heating sources. Eq. (5) tells us that the larger 
grains will tend to dominate both optical extinction and 
submm/mm thermal emission. It is known, however, that 
observations of galactic-scale extinction are weighted to- 
wards widely-distributed, diffuse dust clouds rather than 
localised, dust clumps (Misiriotis & Bianchi 2002). In con- 
trast, larger grains appear to be restricted to compact 
media such as circumstellar shells or pre-stellar cold cores. 
Given the temperature of 10-20 K we measure for the cold 
dust in this study (Sect. 5), we rule out the environs of 
young stars as the primary source of compact submm/mm 
thermal emission within our objects. 

To the best of our knowledge no current astronom- 
ical dataset exists from which we can surmise the pro- 
portion of submm/mm emission emanating from discrete, 
localised sources within the Galactic plane. Such a study 
would allow us to gauge the proportion of submm/mm 
thermal emission arising from compact sources such as 
pre-stellar cores or circumstcllar shells within our targets. 
Utilising COBE data for the Milky Way, Sodroski et al. 
(1997) found that | of 140 and 240 /mi emission comes 
from diffuse HI clouds and is, therefore, unlikely to be 
associated with discrete, compact sources. A comparable 
decomposition of sources at longer wavelengths has not 
been attempted due to the small field of view of current 
submm/mm arrays (2 - 3' for SCUBA and IRAM). This 
is ill-suited to a large-scale survey of the Galactic plane. 

One further consideration is whether grains of the 
same composition are responsible for both optical extinc- 
tion and submm/mm thermal emission. Most grain mod- 
els that have been proposed are composed chiefly of a 
silicate and carbon/graphite component (Whittet 1992). 
The carbon/graphite grains are expected to be more effi- 
cient at absorbing optical radiation (Draine & Lee 1984). 
They will, however, also reach a higher equilibrium grain 
temperature and thus contribute more emission to the 
submm/mm waveband. Temperatures of 15-20 K and 10- 
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15 K, respectively, are expected for carbon and silicate 
grains heated by the general interstellar radiation field 
(Mathis et al 1983). 
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